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Abstract 

We study BPS saturated domain walls in super symmetric SU(N) 
O ■ Yang-Mills theories in the large N limit. We focus on the Seiberg- 

<D Witten regime of N = 2 theory perturbed by a small mass (m < 

0(A/A^ 4 )) and determine the wall profile by numerically minimizing 
its energy density. Similar to the SU{2) wall studied in a previous 
^ . work, the SU(N) wall has a five layer structure, with two confine- 

ment phases on the outside, a Coulomb phase in the middle and two 
transition regions. 
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1 Introduction 



The phenomena of Domain walls, namely, classical configurations that inter- 
polate between degenerate discrete vacua arise in many areas of physics. In 
fact even in high energy physics domain walls show up in very distinct sce- 
narios like in supergravity, in grand unification and in supersymmetric gauge 
theories. In this note we focus on the latter framework. Domain walls in this 
context have been a subject of a rather intensive study in recent years 
- j|l 5|| . Theories with spontaneously broken discrete symmetries are natural 
setups for domain walls. A well known prototype of such a theory is the 
M = 1 Supersymmetric SU(N) Yang-Mills theory. It has a non-anomalous 
Z 2 n chiral symmetry which is spontaneously broken down to the Z 2 group 
by the expectation value of the gaugino bilinear (tr X a X a ). The SYM theory 
thus has N c degenerate discrete vacua, each characterized by a different value 
of the chiral gaugino condensate. 

In supersymmetric theories there is a special class of domain walls, the 
so-called BPS-saturated domain walls which preserve half of the supersym- 
metries of the theory. In terms of the field configuration of the domain wall, 
preserving supersymmetries imposes first-order differential equations of the 
general form 

(Qa ~ ie^Q*) |wall) = 0, (1) 

where the wall is in the xi,x% plane with the fields varying in the normal 
direction x 3 , Q a are the supersymmetry generators, and ip = aig(AW) is the 
phase of the superpotential difference AW between the two vacua separated 
by the wall. The solutions to these BPS equations automatically solve the 
second order field equations of motion. 

Like other BPS-saturated states, the BPS domain walls are more tractable 
and one may reasonably hope for some exact results for such walls even in 
a context of a confining strongly interacting theory. Indeed, the tension i.e. 
energy per unit area of a BPS domain wall is exactly determined by the 
difference between the superpotential values in the two vacua connected by 
the wall without having to find the actual wall profile, 

e=^*=2\AW\ (2) 
Area 

In the M = 1 SYM theory, the superpotential — which acts as a central 
charge for domain walls — is related by the chiral anomaly to the gaugino 
condensate, so a BPS domain wall has tension|| e = ^% |A (tr A Q A a )|. 
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There are two basic ingredients that determine the BPS equation (|T]), 
the superpotential and the Kahler potential. Whereas the exact form of the 
effective superpotential for the M = 1 SYM theory is known, the Kahler 
potential is not fully determined. The effective superpotential is constrained 



by the requirements of holomorphy and flavor symmetry [|16|, [T7|, |18| . Un- 
fortunately, no such constraints apply to the effective Kahler function of the 
theory which controls the kinetic energies of the fields. Prior to our previous 



paper |L9|], all the investigations of this issue have assumed specific Kahler 



functions, only to find that the answer depends on their assumptions. In 
fact, the singularities of a particular Kahler metric led the authors of 0, [3], £|] 
to claim that the SYM theory has an additional chiral-invariant vacuum - 
despite overwhelming evidence that it does not. 

Such problems of determining the BPS equations can be avoided in the- 
ories with more supersymmetries. Indeed, the situation is under much bet- 
ter control for the M = 2 SQCD where the Kahler metric follows from a 
holomorphic pre-potential and the entire low-energy effective Lagrangian is 
completely determined by the Seiberg-Witten theory PP[ . In the Af = 1 
terms, the M = 2 SQCD has an extra chiral superfield in the adjoint repre- 
sentation of the gauge group. Giving this superfield a mass m^O breaks the 
supersymmetry down to M = 1. In the m — > oo limit, the adjoint superfield 
decouples from the low-energy physics and one is left with an effective M = 1 
SQCD. We have therefore decided to study the BPS-saturated domain walls 
in the M = 2 SQCD perturbed by the adjoint mass m. 

The analysis of the the BPS domain walls in the SU(2) Seiberg-Witten 
(SW) theory perturbed by a small adjoint mass, m < A, was performed in 
||T9fl . In the small mass regime, it was found that the SW domain wall has a 
sandwich-like five-layer structure. In each of the two outer layers, the fields 
asymptote to their respective vacuum values. This behavior corresponds 
to the two confining phases of the SW theory characterized by respectively 
magnetic monopole or dyon condensates. In the middle layer, the theory 
is in its Coulomb phase, and the modulus field slowly interpolates between 
its stable- vacuum values; for mass ^ 0, the Coulomb phase is thermody- 
namically unstable in bulk but exists in a layer of finite thickness inside the 
domain wall. The two remaining layers contain transition regions between 
the Coulomb and the appropriate confining phases. This profile was deter- 
mined by solving numerically the BPS equations for each of the wall's layers. 
It was found that the mass is restricted to the region m < A/400. Beyond 
this limit, the transition regions take over the Coulomb phase region and 
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overlap each other. Also, the wall becomes too thin to be analyzed in terms 
of a low-energy i.e long-distance effective theory such as Seiberg-Witten; 
Nevertheless, it was argued that the BPS-saturated domain wall exists for 
any m, small or large. In section 2 we summarize the analysis and results of 

The original SU (2) Seiberg-Witten theory was extended |24], |25| to SU (N) 
gauge groups. Since these theories are more complicated than the 577(2) the- 
ory, so is the analysis of their domain wall configurations. However, one may 
anticipate simplifications by using the large N limit. Indeed such simplifica- 
tions were found in ffl where the large N limit of the SU(N) SW theories 
were addressed. In general the SU (N) Af = 2 SW theory has a moduli space 
of N — 1 complex dimensions. There are N points at which N —1 monopoles 
become massless, which become the vacuum points of the theory when per- 
turbed by a superpotential breaking the number of supersymmetries down 
to one. 

Douglas and Shenker [0] introduced a trajectory along the N = 2 moduli 
space that connects the extreme semiclassical regime to the singular point. 
They showed that along this trajectory the one loop expression for od is exact 
apart from a vanishingly small region around the singularity. To analyze 
the domain wall profiles we introduce a trajectory that interpolates between 
adjacent vacua by analytically continuing the trajectory of KJ. 

The object of this work is to solve the equations and extract the profiles 
of the BPS saturated domain-walls in the SW theory with SU(N) gauge 
group in the large N limit, perturbed by a mass term. The equations for a 
BPS wall appear to be too complicated to solve directly. Instead the wall 
is constructed in two steps: firstly, by regarding only a subset of the BPS 
equations and disregarding the rest, a first approximation of a domain wall 
profile was calculated. Then by a numerical method the configuration was 
iteratively deformed into a state of minimal energy. This two step procedure 
is described is section 3. We found that the final wall configuration is again (i) 
limited to the small values of m and (ii) characterized by a five-layer profile. 
In section 4 we summarize the results and discuss some open questions. 



2 The SU(2) domain wall 

Consider the Seiberg-Witten theory with gauge group SU(2). When per- 
turbed by a superpotential breaking supersymmetry down from Af = 1 to 
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J\f = 2, the moduli space of vacua collapses into two distinct vacuum points. 
At one of these vacua, a massless monopole field condenses, and at the other 
a dyon field condenses. The low energy effective action is thus a different the- 
ory in each of these vacua, and when a domain wall stretches between them, 
different theories rule its different regions. The different theories cannot co- 
exist in the same region of space, and must be spatially separated in order 
for the configuration to make sense. Along the domain wall, three distinct 
phases can be identified: an electric-confinement phase where monopoles ex- 
ist, an oblique-confinement phase with dyons, and a Coulomb phase, away 
from the two vacua where neither monopoles nor dyons are allowed. 

The expectation value of u = tr(0 2 ) serves as the global coordinate of 
the M = 2 moduli space, for the construction of the domain wall. The 
perturbation term for the superpotential is m ■ u. The Kahler metric for 
this coordinate, g ua , is found by numerically calculating the Seiberg-Witten 



elliptic integrals p0| . The resulting metric is a function of u in the domain 
[—A, A], which diverges logarithmically towards both ends, and is rather flat 
otherwise. It was found in ]19[ to be quite accurately approximated by the 
function 

1 i 64A 4 
^H-T^logXT— r ( 3 ) 

The metric for the monopoles is taken canonically to be g M u = g^^ = 1 
and similarly for the dyon fields. For the superpotential, one must look 
separately at different regions of u, and determine it separately for the three 
different phases of the theory. The perturbed superpotential used for the 
electric-confinement phase is 

A 2 - u 

W = m-uA ^MM. (4) 

At the oblique-confinement region a similar superpotential can be written for 
the massless dyons: 

A 2 + u 2 

W = m-u+ — -^ADD. (5) 
— iV2 

Away from the vacuum points, in the Coulomb phase, the only term in the 
superpotential is m • u. This makes sense if the expectation value of the 
monopoles vanishes as we get away from u = —A, and that of the dyons 
vanishes away from u — A. In the construction to follow this condition is 
met, and imposes a limitation on the value of the mass m. 
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Generally, the field profile of a BPS-saturated domain wall is governed 
by eqs. ([]]). In an effective field theory without higher-derivative terms in its 
Lagrangian, the explicit form of these equations is 

^=e-„«— (6) 
dx z 9 dA> { ' 

where A l (x$) runs over all the scalar fields of the theory. For the theory 
at hand, this means u, M, M, D and D. Fortunately, in the Seiberg- 
Witten regime of the theory (small m), the monopole and the dyon fields are 
segregated to different regions of space, so for any particular stretch of £3 we 
may eliminate some of the fields from the BPS eqs. (^). 

Altogether, the wall has five distinct layers: On the left side, the u, M and 
M fields asymptote to their vacuum values in the electric confinement phase 
and there are no dyons. Next, there is a transition region where monopole 
fields turn off. The third, middle layer is in the Coulomb phase, where both 
monopoles and dyons are absent and u varies slowly from A 2 to —A 2 . Next, 
there is another transition layer where the dyon fields turn on. Finally, on 
the right side, u, D and D fields approach their vacuum expectation values 
in the oblique confinement phase. 

Consider the middle layer first. The only important field here is u, hence 
the BPS eqs. reduce to 

= m *e^g m . (7) 

CIX3 

Without loss of generality, we assume real positive m and A, hence e l<p = — 1 
and real u(xs). For the approximate metric (|3]), the exact solution is 

2 , A 2 + u u 64e 2 A 4 

16vr mx 3 = log^-^ - ^log^^ (8) 

and the do main- wall has finite thickness (as defined in |19[]) of about w = 
O.OGOSm^ 1 . Working numerically with the exact metric yields a slight ad- 
justment to this value, w = 0.0625m _1 . 

The two transition regions are isomorphic, so it is enough to consider the 
transition between the electric-confinement and the Coulomb phases. The 
relevant fields here are u, M and M, the superpotential is @ and thus the 
BPS equations are 
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Figure 1: Field profiles for the Seiberg-Witten domain wall, (m = A/2000). 
(Figure 5 of @) 

dM A 2 -u* ~ 

M , (9) 



dx 3 i\/2A 
dM A 2 - u* 



-M*. 



dx 3 iy/2A 

As in the Coulomb phase, the u field takes real value while the magnetic 
gauge symmetry allows us to specify real M and imaginary M = iM. The 
boundary conditions for eqs. ([)]) amount to requirements that 1. all fields 
must asymptotically approach their vacuum expectation values for x —>■ — oo, 
and 2. the monopole fields M, M must disappear towards the middle of the 
wall so that the solution can be matched with the previous solution of the 
Coulomb phase. Of course, analytically the monopole fields cannot vanish 
exactly, but numerically they do become negligibly small at x 3 = — drop 
to less than 1% of their vacuum values — provided m is smaller than about 
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A/400. 

The numerical wall profile for m = A/2000 is shown in figure ([I]); note the 
transition regions occupying a considerable fraction of the wall width despite 
rather small m, although this fraction becomes smaller for even smaller values 
of the perturbing mass parameter. On the other hand, for larger perturbing 
masses, the transition regions spread towards the mid-wall point and start 
overlapping each other. Both the monopole and the dyon fields appear to 
be present in such an overlap, a situation quite untractable in terms of an 
effective field theory. In fact, this is just a symptom of a more general 
problem: A larger m makes for a thinner wall, which can never be adequately 
described in terms of a low-energy — i.e., long-distance — effective theory 
theory such as Seiberg-Witten. 



3 The large- N domain wall 

The A/" ~ 2 limit of the SU(N) SYM is generally similar to the Seiberg- 
Witten N = 2 theory, but instead of a single complex modulus u = tr(0 2 ) 
we now have N — 1 moduli. The eigenvalues 1 , . . . , <fr N of the scalar field <ft 
(viewed as an A x A matrix) provide a redundant coordinate system for this 
moduli space (note J2j = 0). At A" discrete points in this space [I], 

(f) = e™ T ' N x 2Acos ^ ^ 1/2) , r = 0,l,...,(JV-l), (10) 

the maximal number of magnetic monopoles or dyons become massless. Such 
massless monopoles/dyons acquire non-zero vacuum expectation values when 
the Af = 2 theory is perturbed by the superpotential 

^=yEW 2 , (11) 

3 

(the M = 1 mass term) and the theory enters a confining phase (electric or 
oblique). All vacua other than (|10D become unstable. 

Following M. Douglas and S. Shenker Jjj], let us Fourier transform the 
eigenvalues (fP according to 

t = A|; /C08 ^_VH>, (12) 
¥ = " E V cos £Mz_M. (13) 

k=l 
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In terms of the new moduli t 1 , . . . , t^ -1 , the N stable vacua are located at 
t 1 = 2he mT l N and t 2 = t 3 = ■ ■ ■ = t k = 0. Let us focus on the r = vacuum 
(electric confinement) where the light fields comprise the t k moduli and the 
magnetic monopoles Mi, Mi (I = 1, . . . , N — 1). The superpotential for these 
fields combines the perturbative term (pi]) and the monopole mass terms 
YjO-fi^MiMi. Expanding the latter around t k « 2A5j. i1 and taking the large 
iV limit [||], we arrive at 

W = ^E(^) 2 + £sin^(t fe -2A4,i)M,M,. (14) 
4 k k,i ^ 

Among other things, this superpotential governs the monopole condensa- 
tions; solving for dW = 0, we find 

MM) = -2AAmsin^. (15) 

3.1 Domain Wall: First Approximation 

We would like to construct a BPS-saturated domain wall between two adja- 
cent vacua (pTDj) , say r = and r = 1. Both vacua have t 2 = ■ ■ ■ = t N ~ l = 0, 
so as a first approximation to the domain wall's profile, we assume all the 
t h>1 moduli to maintain zero values throughout the wall while the t 1 mod- 
ulus evolves from t 1 = 2 to t 1 = 2e m ^ N . (Henceforth we shall use A = 1 
units). This 'trajectory' through the moduli space is related to the 'scaling 
trajectory' 



^( g ) = e s / Jv x2cos 7T(j 1/2) (16) 

of ref. pj via analytic continuation to imaginary s that runs from to iri 
(t 1 = 2e s / N ). Along this trajectory, the moduli metric is diagonal in the t k 
basis (in the large N limit only!), 

Re(FG*) 

9kk> — d k,k f x q : ' i 1 7 ) 

S7T sin ^ 

where F(s, j^) and G(s, ^) are defined in eqs. (5,9-10) of |J; for k <C N, 
F w G ~ 1 regardless of s. For larger k = O(N), F and G — and hence 
the metric g^ k — become s dependent, but fortunately, the detailed nature 
of this dependence is not germane for the present analysis. 
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Similarly to the SU (2) case, we expect the domain wall to have the five- 
layer structure with two well-separated transition regions for sufficiently small 
perturbing mass m. Consequently, there should be no monopole fields for 
13 > or dyon fields for x 3 < 0. Physically, the wall is left-right symmetric, 
so let us focus on its left half and thus dispense with the dyon fields. The 
remaining moduli and monopole fields are governed by the BPS equations 

*£ ■ *S • 

—-L - e^^L, (19) 



dx 3 dMi 
dM,* dW 



dx 3 dM t 



(20) 



In the Coulomb phase domain where the monopole fields have very low values, 
the right hand sides of these equations vanish, except for the the first eq. fll8|) 
for the t l modulus. Thanks to the diagonality of the moduli metric ([T7|), this 
means that t 1 is the only non-constant field in the Coulomb domain, in perfect 
consistency with the t k>1 = trajectory assumption. Unfortunately, in the 
middle of the transition region where the monopole fields neither vanish 



nor take vacuum values, the right hand sides of eqs. (|18 ) for the moduli 
t 2 , t 3 , . . . , t^™ 1 do not vanish, which means that the t k>1 = wall is not 
exactly BPS. 

Nevertheless, as our first approximation, we would like to freeze all t k>1 = 
and solve the BPS equations for the ^(x^) and all the monopoles Mi(x 3 ) 
and Mi(xs) analytically. In the following subsection, we numerically calculate 
the deviation of the true BPS wall from this trajectory; we shall see that 
such deviations are negligible in the Coulomb domain but significant in the 
transition region, especially its outer edge. For the overall profile of the wall 
however, t k>1 w is a good approximation. 

For the two vacua r = and r = 1, AW = iV 2 m(e 27n — 1) ~ 2irNmi, so 
let us re-phase W — > —iW and make the eqs. (|18D [for the T 1 ] and (|ITH2~0l) 
real. The boundary conditions are t l = 2 and eqs. flT5p for the monopole 
fields at x 3 = -00, while for x 3 = 0, t 1 = 2e ni/2N w 2 + iri/N and the 
monopoles must vanish. The BPS solution can be generally written as 



Mi = Mi = iJ2Nmsm(^) x exp(-/i sin ^) , (21) 
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1 2is 1 . , 

t 1 = 2 + — , (22) 

with s x (x3) and h(x 3 ) developing according to differential equations 

<)S> 2n 2 Nm[N -2 J2^ 2 (-)e- 2hsin % (23) 



dx 3 1 i N 

dh 2s 1 
dx~ 3 ~ ~W 



(24) 



Physically, s 1 simply keeps track of the coordinate t 1 of the Af = 2 moduli 
space, while h governs the decay of the monopole condensates away from the 
confining vacuum domain; note different monopole condensates Mi decaying 
at different rates thanks to the sin(^) factor in the exponent in eq. ([2"T|). 

The boundary conditions are s 1 — > and h — > at — oo, and s 1 = n/2 
a.t x 3 = 0. Also, all the monopole condensates must become vanishingly 
small towards £3 = 0, thus h{ji/N) ^> 1 at X3 = 0. Physically, the last 
requirement imposes a lower limit on the domain wall's width and hence an 
upper limit on the perturbing mass m. Indeed, in light of eq. fl24|) we estimate 
h(x 3 = 0) ~ x f ' hence the wall should be much wide than wo ~ 4N 2 /ti 2 . 
On the other hand, eq. ( ^3] ) lets us estimate 2 r n 2 N 2 m x -f ~ f , hence 

w zttt- ■ (25 

Consequently, to maintain spatial separation between the monopole and the 
dyon fields, we need 

m < (26) 
8N 4 K 1 

(in A = 1 units). 

Given small enough m, equations (^3]- |2^) can be easily solved numerically. 
Figure (^) displays the resulting field profiles through the transition region of 
the wall, as calculated for A^ = 51, m = 2.2 ■ 10~ 9 . Note different monopole 
condensates decreasing at different rates but all reaching negligible values in 
the Coulomb domain (the right side of the figure). 



3.2 Correcting the Trajectory: Minimizing the Energy 

The net energy density of a BPS-saturated domain wall comprises potential 
and kinetic terms, each contributing \AW\ to the total. The approximately 
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Figure 2: First approximation to the transition region of the domain wall, 
between the electric confinement phase on the left and the Coulomb phase 
on the right. One can see how the monopole condensates disappear towards 
the middle of the wall (the right side of the picture), while the s 1 coordinate 
moves from up to | at the center (far right, out of the frame). The hor- 
izontal axis is x% in units of 1/A, shifted so that the zero point is deep in 
the confining phase on the left. The scale for the monopole condensates is 
shown on the left; the scale of s 1 is 2500 times large, thus s 1 = 0.75 at the 
right edge of the frame. 
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but not exactly BPS domain wall of the previous section has correct kinetic 
energy, 

9i] d 3 A%A3 = gud^dst* 1 +J2(d 3 M l d 3 M; + d 3 M l d 3 M' l 

i 

= dst 1 —- + Y d 3 Mi-— + d 3 M t ^ 

dt 1 ^ 3 l dM l l dMj 



thus 



dx 3 d 3 W 



AW. 



(27) 



Its potential energy however is higher than BPS because of un-satisfied BPS 
equations fll8 ) for the t k>1 . Indeed, 



-pot 



dx 3 g lJ 



dW dW* 



dA OA, 



dx 3 
Jdx 3 



Y,9 

k 



kk 



dW dW* 



dt k dt 
1 dW 



F + Y {dMidMf + dMidM*' 



kk 



dW 



dt k 



aw + dx 3 y g 

J k^l 



kk 



dW 



dt k 



(28) 



Altogether, the last term in eq. ( f28|) is the excess energy density over the BPS 
domain wall due to violation of some of the BPS equations. Spacewise, this 
excess energy comes solely from the transition region (dW/dt k = in both 
Confining and Coulomb phases), thus our approximate solution is accurate 
for the Coulomb phase region — which constitutes most of the wall width 
and provides most of its energy in the small m limit. However, the transition 
regions of the wall are rather interesting even when they are narrow, so we 
would like a better approximation to the field profiles. In particular, we need 
to get rid of the wrong assumption about t k>1 and calculate their actual 
profiles. 

To obtain such a better approximation, we minimize the wall's net en- 
ergy density as a function of its field profile. Basically, we first discretize x 3 
and analytically evaluate the variational derivatives of the energy with re- 
spect to all possible lattice field variations, then numerically evolve the field 
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configuration in the direction of the steepest energy descent. The boundary 
values are fixed to vacuum values at the left boundary and all fields except 
t l are zero at the right boundary. We start with the approximate wall profile 
obtained in the previous subsection, then evolve the lattice fields until we 
reach the (numeric) minimum of the wall's energy. 

Our calculations are somewhat eased by the symmetries of the superpo- 
tential (|l4j) and hence of the exact solution to the BPS equations, 

M, Mi M N i = M N , and t k = Veven k. (29) 

Naturally, we hardcoded these symmetries into our minimization procedure. 
Also, for convenience of taking the large N limit, we use rescaled moduli 
fields 

s k = ^(t k -26 kA ); (30) 

note consistency with the s 1 variable used in eq. (p2[). In terms of the inde- 
pendent Mi{x^) and s k {x^) ) the energy density we need to minimize is 

dx 3 p (31) 



where 

PM = E {—^(^ + — E Mfsm— ~ -—8 k A 

oddfc [ 4 9kk \ i=i ^ 1 I 

+ (7V E /2 UdM f + H M f ( |> fc sin (32) 



and 



9kk ~ o — : — W (33) 
87r sm 



since Re(FG*) « 1 in eq. flT7|).g 



2 To be precise, this approximmation works very well for k <C -/V but for A: = O(N) one 
should in principle use a more accurate approximmation 

m«n - (i - jl) (i - 1) (i - £) + \ £ x i„ g i . (3 4, 

Fortunately, numerically this expression is close enough to 1 in the region of interest to 
justify the approximmate metric (p3|). 
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The technical details of the minimization procedure being rather boring, 
let us simply present our results. Figure (|3|) depicts the monopole conden- 
sates profiles — which have not changed dramatically during the energy- 
minimization procedure, cf. figure (|J), except their transition regions be- 
came a bit narrower. This is easily understood in terms of BPS eqs. ([i8Hl9|): 
the additional terms involving t k>1 moduli make for steeper decrease of the 
monopole fields. Figure (^) displays the moduli profiles s 1 and also s 3 , s 5 , 
etc. As expected, s 3 , s 5 , . . . vanish in the Coulomb phase on the right, but 
not in the transition region between the confining and the Coulomb phases. 
For example, at its peak in the middle of the transition region, the s 3 reaches 
about 18% of the value of the s 1 at the same point. For other odd k, the 
s k reach their peaks at about the same point, reaching a value roughly pro- 
portional to 1/k. The reason for such behavior is the superpotential fll4|): 
for large k, sm(iikl/N) oscillates faster between positive and negative values, 
leaving an average of order 1/k. 

The width of the transition region can now be measured by looking at 
either monopole condensates or moduli s h . The highest peak of the s 3 field 
- where its value is above 50% of its maximum — corresponds to the region 
where the highest monopole mode M^-i)/2 decreases from 85% of its vacuum 
value down to 15%. Other monopoles are wider and recede less steeply. As 
for the higher moduli fields, the s fc>3 have wider peaks, for example twice 
as wide for s 21 than for s 3 and slightly more than that for yet larger k. We 
see no obvious way to account for this fact, other than to say that it is a 
numerical result. In any case, the high modes are not that much wide wider: 
merely by a factor of 2 rather than N or even v^A- 

Interestingly, on the immediate right of the transition region, the s h>1 
moduli drop below zero and develop small negative values before they finally 
disappear in the Coulomb region. The origin of this behavior is unclear, but 
this is definitely not an artefact of the boundary condition. This is verified 
by pushing the right boundary of the lattice region further right. 

Finally, we would like to mention that when all the s h are very small, 
they have comparable magnitudes, for example s 3 ~ 0.65s 1 . In other words, 
at the left edge of the wall, where the moduli fields just begin to shift away 
from their values in the confining vacuum, they move in a direction at rather 
wide angle to the s 1 "axis" characteristic of the Coulomb phase in the middle 
of the wall. This rather unexpected geometry of the wall's 'trajectory' of the 
moduli space indicates that the present analysis does not extrapolate well 
into the regime of not-too-small (A/" = 2)-breaking mass m for which various 
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Figure 3: Field profiles for the monopoles at the final domain wall configura- 
tion, N — 51 and m — 2.2 • 1CT 9 . The highest and steepest line is M 25 , and 
the lowest M±. 
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Figure 4: Field profiles for the moduli s k (odd k). The s 1 will continue to 
rise throughout the Coulomb but the plot is cut at x — 3200 to avoid lines 
sticking out of the frame. The tallest line after s 1 is s 3 , then s 5 , etc. 
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regions of the domain wall start overlapping each other. 
3.3 Conclusion 

The structure of the domain wall worked out in this work is in many ways 
reminiscent of the SU (2) domain wall. Both exhibit a similar five-layer pro- 
file, and both depend strongly on the smallness of the perturbation mass pa- 
rameter to. The devil is in the details and the SU (N) theories exhibit many 
new features associated with multiple monopole/dyon condensates and mul- 
tiple moduli fields. There are many more BPS equations now, all entangled 
with each other, so calculating the field profiles of the wall calls for numerical 
lattice calculations. Even the BPS route that the moduli fields take on their 
way from one stable vacua to another is highly turns out to be non-trivial. 
In the SU (N) model, the route between two vacua in the moduli space is not 
trivial, and distinct functions describe the evolutions of different monopole 
fields. 

The large N approximation helps us to get explicit expressions for the 
Kahler metric and for the superpotential, but unfortunately, at the end of 
the day we are again faced with the same limitation as in the 577(2) domain 
wall, namely that the wall is valid only for a small perturbation parameter 
to. In fact, the limitation here is more severe, as to must be smaller than 
0(N~ 4 ). For such small to, the five-layer picture is adequate. The work 
shows that the dominant field in the Coulomb phase is t 1 , and the other 
fields either asymptote to zero, or remain very small. It is not possible, 
however, to disregard the other t's in the transition region. As to gets larger 
and the two transition regions overlap, on top of the difficulty from having 
coexisting monopole fields from two different vacua, we cannot even say along 
what route in the moduli space of the M = 2 theory the domain wall lies. 
The wall derived here, however, is a legitimate domain wall, and it may be 
interesting to study its low energy excitations in the future. 
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